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ability  of  immediately  meeting  customer  requisitions  on  high  priority  items; 
and  c)  meeting  previously  unfilled  demands  on  high  priority  items  as  quickly 
as  possible. 


We  develop  a  chance  constrained  program  for  this  problem  that  will 
•yield  decision  rules  expressing  the  decision  variables  (quantities  of  each  item 
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to  order  at  each  period)  that  are  consistent  with  the  constraints  above  and 
which  maximize  the  sum  of  the  probabilities  of  meeting  customer  requisitions 
on  all  items  over  a  specified  number  of  periods.  We  then  develop  an  equivalent 
deterministic  program  whose  solution  will  determine  the  coefficients  of  these 
decision  rules. 
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Chapter  1 


Introduction 


1.1  Overview 

We  will  address  inventory  control  for  consumable  items  (verses  re¬ 
pairable  items)  at  Air  Force  base-level  supply  stores.  Our  objective  is  to  manage 
the  inventory  at  the  supply  store  so  that  we  meet  as  many  customer  requisitions 
as  possible  for  n  items  over  N  periods  subject  to  a  budget  constraint.  Customer 
requisitions  and  re-supply  leadtimes  are  stochastic  random  variables  considered 
stochastically  independent.  Further,  we  want  to  be  confident  of  meeting  cus¬ 
tomer  requisitions  on  high  priority  items.  We  assume  that  requisitions  we  can’t 
meet  from  stock  on-hand  can  be  supplied  later  (i.e.  complete  backlogging),  but 
we  want  to  guarantee  that  unfilled  requisitions  for  high  priority  items  will  be 
filled  as  soon  as  possible. 

In  this  report  I’ll  start  with  a  single-period  model  that  isn’t  very 
realistic  but  whose  structure  allows  solutions  to  be  developed  more  easily,  then 
move  to  a  more  complicated  model  which  is  much  closer  to  reality.  The  first 
model  will  seek  to  minimize  Expected  Backorder  Costs  times  a  weighting  factor 
set  for  each  item. 

The  second  model  will  seek  to  maximize  the  sum  of  the  probabilities 
of  meeting  customer  requisitions  on  n  items  across  N  periods  subject  to  chance 
'  ’constraints  concerning:  a)  dollars  obligated;  b)  meeting  requisition  demands 
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for  important  items;  and  c)  meeting  unfilled  requisition  on  important  items  as 
quickly  as  possible. 

1.2  Background  on  the  Current  Air  Force  System 

Currently  the  Air  Force  bases  use  a  software  package  called  the  Stan¬ 
dard  Base  Supply  System  (SBSS)  to  compute  the  leveling  parameters  for  each 
item  in  their  consumable  inventories.  The  SBSS  sets  these  parameters  for  each 
item: 


r,  Reorder  Point-when  the  inventory  position  of  item  i  is  below  r,  then  an 
order  is  placed 

<7,  Quantity  to  Order-any  orders  placed  will  be  for  this  quantity 

Base-level  purchases  for  consumable  items  are  funded  through  a 
“stock  fund.”  New  orders  can  be  placed  to  replenish  the  inventory  only  if 
enough  funds  are  available  in  the  “stock  fund.”  When  a  customer  “buys  '  an 
item  from  the  base’s  inventory,  money  for  the  item  is  transferred  from  one  of 
his  funds  (usually  O&M-Operations  and  Maintenance)  to  the  “stock  fund.” 

While  the  formulas  used  to  set  these  parameters  in  the  SBSS  are 
designed  to  reduce  long-run  costs  for  the  inventory,  they  ignore  any  constraints 
on  the  solution.  The  parameters  are  computed  independently  for  each  item, 
based  only  on  the  items  past  demand  pattern,  cost,  and  estimates  of  ordering 
and  holding  costs.  These  parameters  are  intended  to  be  used  in  a  continuous 
review  inventory  system,  i.e.  the  inventory  level  of  each  item  i  is  continuously 
monitored  and  when  the  level  drops  below  r,,  an  order  for  <7*  units  of  item  i 
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should  be  placed.  But,  since  these  orders  must  be  funded  from  the  “stock  fund,’" 
often  its  impossible  to  place  this  order.  Consequently  the  parameters  the  SBSS 
computes  are  often  unusable  when  daily  decisions  are  made  on  what  quantities 
of  which  items  need  to  be  order  while  staying  within  budget  constraints. 


Chapter  2 


A  Single  Period  Model 


In  this  chapter  we  will  consider  a  formulation  of  the  Inventory  Control 
problem  for  a  single  period  of  length  L  days.  The  objective  will  be  to  determine 
a  single  order-up-to  parameter,  yi  fbr  each  item,  so  that  if  the  starting  inventory 
is  yi  for  each  item,  then  the  sum  over  all  items  of  Weighted  Expected  Backorders 
will  be  minimized  for  the  period.  These  y,  must  be  chosen  so  that  the  value  of 
the  inventory  does  not  exceed  a  budget  constraint  of  7maxdollars. 

2.1  Notation 

L  Length  of  Period  in  days 

0  Random  Variable  for  the  number  of  customer  demands  for  the  period. 
[UNITS] 

fi(0)  Probability  density  function  for  the  number  of  demands  of  item  i  for  the 
period 

Ci  Unit  Cost  of  item  i  [per  UNIT] 

Imax  Budget  Constraint  in  $’s  for  on-hand  inventory 
y<  VARIABLE-optimal  starting  on-hand  inventory  for  item  i.  [UNITS] 

7r,  Weighting  Factor  for  each  backorder  of  item  i  [per  UNIT] 
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Ei{yi )  Expected  number  of  backorders  for  item  i  for  the  period  if  the  starting 
inventory  is  j u  [UNITS] 

2.2  Assumptions 

•  Single  Period  model  of  length  L  days 

•  Stock  is  raised  to  t/;  before  beginning  of  period 

•  No  ordering  of  stock  is  permitted  during  period 

•  Maximum  investment  in  inventory  is  /maIdollars 

•  Demand  for  each  item  is  stochastically  independent  of  the  other  items' 
demands 

•  The  distribution  of  demands  for  each  item  i  can  be  accurately  approxi¬ 
mated  by  a  known  finite  discrete  distribution. 

This  single-period  model  is  formulated: 

min  £•=!*.■£.(!/.)  (2.1) 

S.t.  £"=1c,t/t  ^  v 

yi  >  0  Vi  =  1 , . . . ,  n 

yi  is  an  integer  Vi  =  1, . . . ,  n 

with  Ei(yi)  =  the  expected  value  of  non-negative  (6  —  yi) 


Ei(yi)  =  myi(0-yt)fi(0)de 
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2.3  A  Continuous  Relaxation 

If  we  drop  the  constraint  that  the  computed  y<  be  integers  and  instead 
treat  them  as  continuous,  our  problem  becomes 


n 

min  S>.£,(s.) 

«=1 

s.t.  E”=ic«y,  <  v 

Vi  >  0  Vi  = 


(2.2) 


The  expected  number  of  backorders  of  item  i  for  the  period,  given  y,, 
denoted  Ei{yi)  is 


Eiivi)  =  -  yi)M0)de.  (2.3) 

Note  that  £,is  truncated  for  0  values  that  can  be  met  by  on-hand 
inventory,  i.e.  for  0  <  y,. 

By  replacing  ( 0  —  y, )  by  max{0,  ( 0  —  y*)},  denoted  [0  —  y,]+,  we  can 
rewrite  E{  integrating  over  all  nonnegative  values  of  0. 

Ei  =  f°°[0  -  yi]+ fi{0)d0  (2.4) 

Jo 

Let  pi'O,  =  Probability{a  units  of  item  i  are  demanded}.  By  the 
assumption  that  the  demand  can  be  approximated  by  a  finite  discrete  distri¬ 
bution,  there  exists  fc,  for  each  i  such  that,  =  1,  then 


Ei(Vi)  =  Earl  (Pi,a  *  [0a  ~  2/,]  +  ) 


So  our  model  can  be  formulated 
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min  J2  {^^=1  (p<>  *  [e»  “  2/.]+)} 

zr=i  ^  <  v  (2-5) 

y,  >  o  V*  =  l,...,n 

With  the  E{  expression  written  in  this  form,  the  problem  is  still  a 
non-linear  problem  (since  the  max  function  is  non-linear).  However,  since  [>]+ 
is  convex,  the  objective  function  is  convex. 

Now,  since  the  indices  of  the  objective  function  are  now  constant, 
i.e.  not  dependent  on  the  this  model  can  be  coded  for  optimization  by  any 
non-linear  optimization  package  that  allows  non-smooth  objective  functions. 

However,  we  can  get  better  results  if  we  use  the  following  completely 
equivalent  Linear  Programming  Model. 
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2.4  An  Equivalent  LP  formulation 

The  only  non-linear  term  in  the  objective  function  is  \0  —  y,]+.  This 
maximum  function  can  rewritten 

max{O,0  —  yi)  =  |  {\e  -  y,|  +  (0  -  Vi)}  (2.6) 

Now,  |#  —  y{\  is  the  only  nonlinear  term.  To  replace  by  a  linear  term. 

let  0a  -  y,  -  £fo  +  6~a  =  0  (2.7) 

where  S  s  are  non-negative  slack  variables  and  0a  represents  demands  of  a  units. 

Lemma  1  y'  for  the  following  formulation  will  also  be  optimal  for  (2.5). 

k, 

£ 

Qf=  1 

Proof:  For  each  z,  a  combination,  if  pt>  >  0,  then  at  the  min  of  this  function, 

'  ^i,a  =  0.  i.e.  at  least  one  of  the  6' s  is  zero.  Otherwise  3c  >  0  such  that 

Ka  ~  Ko  ~  c  and 
Ko  =  Ka~€ 

satisfy  condition  2.7  and  produce  a  lower  minimum  value  for  the  objective 
function  2.8.  Since  one  of  the  <Fs  are  zero,  \0a  -  y,|  =  6fa  -f  S~a.  So  2.8  is 
equivalent  to  2.5  for  z,  a  combinations  such  that  pj  a  >  0. 

On  the  other  hand,  for  z,  q  combinations  where  pt  o  =  0,  then 


Pi, a  [<5,^  +  6~a  +  0Q  -  yi  | 


Pi<<*  *  [$!"»  +  6i,a  +  0  -  y»]  =  o  =  pua  *  [0  -  y,]+. 
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and  2.8  is  trivially  equivalent  to  2.5.  QSV 

So  we  can  reformulate  our  initial  problem  as  the  following  completely 
linear  programming  problem. 


min 

2 

=1  Pi, 

+ 

+ 

oa- 

-S.]} 

s.t. 

< 

fmai 

Vx 

+  Ka  ' 

= 

0Q 

for  a 

= 

1,. 

..,ki 

for  i 

= 

1,.. 

■ , n 

i,ax  Vi 

> 

0 

for  i 

= 

1,.. 

. ,  n 

Conclusion:  If  we  approximate  our  demand  distribution  with  a  finite 
discrete  distribution,  our  initial  non-linear  formulation  (2.5)  can  be  written  as 
the  completely  equivalent  linear  programming  problem  (2.9).  In  fact,  if  we 
compute  the  solution  with  an  adjacent  extreme  point  algorithm,  then  by  the 
LIEP  theorem,  *at  each  iteration  (basic  solution),  at  most  one  of  the  <^a,<5"Qwill 
be  positive. 


'Prof.  Charnes’  Linear  Independent  Extreme  Points  (LIEP)  Theorem  (1950),  see  [5]. 


Chapter  3 


A  Chance  Constrained  Multiperiod  Model  for 
Inventory  Management 

3.1  Problem  Statement 

The  objective  is  to  maximize  the  sum  of  the  probabilities  of  filling 
all  customer  requisitions  at  a  base  supply  store  for  n  inventory  items  over  the 
specified  N  periods  through  the  determination  for  each  period  of  (stochastic) 
decision  rules  from  the  class  of  “informationally  feasible”  rules  linear  in  the 
requisitions.  The  decision  variables  expressed  by  these  rules  will  be  the  quan¬ 
tity  of  each  item  to  order  for  each  period  which  will  satisfy  probabilistic,  e.g. 
“chance”,  constraints  concerning:  a)  dollars  obligated;  b)  meeting  requisition 
demands  for  important  items;  and  c)  unfulfilled  customer  requisitions  for  im¬ 
portant  items. 

We  will  concentrate  only  on  consumables  items  (as  distinct  from  re¬ 
pairable  items).  The  Air  Force  refers  to  the  items  we  will  concentrate  on  as 
“Budget  Code  1”  items  of  “XB3”  items.  These  items  are  ordered  by  the  base 
supply  store  and  paid  for  from  the  base’s  “stock  fund”  (managed  by  the  base). 

3.2  Assumptions 

•  Customer  requisitions  (called  “Sales”)  to  be  filled  each  period  are  random 
variables  with  known  joint  distribution  functions. 
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•  Orders  in  each  period  are  placed  once  by  the  supply  store. 

•  Time  of  arrival  of  orders  (“Receipts”)  from  time  of  order  is  stochastic,  i.e. 
the  probabilities  of  arrival  in  later  periods  are  known.  When  an  order  for 
an  item  arrives,  the  quantity  is  exactly  what  was  ordered.  Further  the 
stochastic  arrivals  of  orders  are  assumed  to  be  stochastically  independent, 
both  from  other  orders  of  the  same  item  and  orders  of  different  items. 

•  Sales  for  each  period  that  cannot  be  met  from  stock  on-hand  at  the 
beginning  of  the  period  can  be  supplied  later. 

•  We  can  identify  a  subset  of  important  items  for  which  we  want  to  insure 
with  high  probability  meeting  customer  requisitions. 

•  The  cost  of  each  item  is  constant  throughout  the  N  period  horizon. 

•  Money  for  each  order  is  obligated  in  the  same  period  in  which  the  order 
is  placed.  Funds  for  the  requisitions  (sales)  in  a  period  are  received  in 
the  period  the  requisition  is  placed,  even  if  the  requisition  is  filled  in  a 
later  period. 

3.3  Formulation 

Let 

c,  [dollars  per  unit]  Cost  of  each  unit  of  item  i. 

Ima x  [dollars]  The  upper  budget  limit  on  the  amount  of  money  that  can  be 
invested  in  the  inventory. 
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S;(t)  [units]  Customer  requisitions  (“Sales”)  of  item  i  during  period  t  (random 
variable) 

Ri(t)  [units]  Receipts  of  item  i  up  through  period  t.  These  units  are  available 
to  satisfy  unfilled  requisitions  through  period  t  —  1  (random  variable) 

Yi(t)  [units]  Quantity  of  item  i  to  order  at  period  t  (Decision  Variable) 

Ii(t)  [units]  Inventory  on-hand  of  item  i  at  beginning  of  period  t.  This  quantity 
is  available  to  satisfy  Sales  for  period  t. 

•  for  i  =  1 , . . . ,  n 

•  for  t  =  1, . . . ,  N 

3.4  Formulation 

We  want  to  maximize  the  sum  over  all  the  N  periods  and  over  all  the 
n  items  of  the  probabilities  of  satisfying  requisitions  (“Sales”),  i.e. 

maxEEW)  tSi{t)) 

t=i  «=i 

subject  to  the  various  constraints  described  generally  above. 

These  constraints  contain  stochastic  elements  both  with  regard  to 
time,  e.g.  as  in  receipt  of  orders  placed,  and  with  regard  to  conformance  of 
ordering  with  Air  Force  stockage  policy. 
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3.5  Chance  Constraints  and  Their  Random  Variables 

Ii(t)  is  a  complicated  function  of  the  previous  orders,  receipts  and 
sales.  We  suppose  that  individual  receipts  are  in  the  same  amounts  as  orders. 
However,  we  assume  that  the  number  of  periods  between  placing  an  order  and 
that  order  being  available  to  fill  customer  requisitions  is  a  discrete  random 
variable  with  known  probabilities.  (These  probabilities  may  be  different  for 
each  item  in  the  inventory.)  If  an  order  Yi(t)  is  placed  its  receipt  cannot  be 
used  to  satisfy  customer  requisitions  (sales)  until  period  <4-1,  even  if  the  order 
is  received  in  period  <.  Of  course,  it  may  not  be  received  until  period  <  +  1, 
<4-2,  etc. 

We  proceed  to  determine  the  stochastic  receipts  of  an  item  at  a  par¬ 
ticular  period  in  terms  of  the  amount  and  period  of  order  placement. 

For  <  >  a,  define 

p,(<  —  5)  =  Prob(Order  ^(.s)  arrives  precisely  in  period  <). 

It  should  be  noted  that  p,(<  —  s)=0  for  <  —  s  >  say  4. 

Let  Ri(Yi(s)  ,<)  be  the  amount  received  by  period  <  of  an  order 
Vi(s)  placed  in  period  s. 

Then  we  set 

0,  s  >  t 

Yi(s)  with  probability  p,(<  —  s)for  s  <  t 
where?  p,-(-)  are  known  probabilities. 


Ri(Yi(s)  ,<) 


< 
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Since  the  random  variables  of  the  arrivals  of  receipts  (of  the  same  or 
different  items)  are  assumed  stochastically  independent,  we  can  write  Ri(t), 
the  cumulative  receipts  of  item  i  through  period  t,  as  the  sum  of  independent 
Ri(Yi(s)  ,t )  random  variables, 

£.(0  =  X>(KOO  ,<) 

4=1 

In  other  words,  since  Ri(Yi(s)  ,  f)  is  either  zero  or  the  amount 
Yi(s)  ordered,  Ri(t)c an  only  take  on  values  that  are  a  im  of  one  or  more 
of  the  previous  ^(-s)  ’s. 

For  example, 


Ri(  i)  = 


^(2)  = 


i?i(3)  = 


Let 

Zi(t)  = 

hi  = 


Value  with  probability 


o,  (l-Pi(l)) 

Vi(l),  ft(l) 

Value 

probability 

o, 

(1  ~Pi(0))(l  -P,(l)) 

Yi(  1), 

P«(l)(l  -  P.(0)) 

m 

P.(0)(1  -p,(l)) 

K(l)  +  K(2), 

p.(l)p«(0) 

Value 

with  probability 

0 

(l-p,(2))(l-p,(l))(l-p, 

15(1), 

p,(2)(l-p,(l))(l-p,(0)) 

>5(2), 

(1  ~  P«(2))p«(l)(l  —  p«(0)) 

15(3), 

(1  —  P*(2))(l  —  P«(l))p<(0) 

y<(l)  +  ^(2), 

P«(2)p»(l)(l  —  P*(0)) 

15(1)  + 15(3), 

P»(2)(l  —  P«(l))Pi(0) 

K(l)  +  K(2), 

(1  -p,(2))p,(l)p,(0) 

tf(l)  + 15(2) +  15(3),  P:(2)p,(l)p,(0) 

(1  -Pi(j)) 

Piti) 

-p.(j)) 


Then  in  general, 
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Ri(t) 


Value 

probbability 

0 

2i(t) 

Y(  1) 

Yi(  2) 

Pi,2Zi{t ) 

*(<) 

Pi,lPi,  2Zi(t) 

Yi(  1)  +  tf(3) 

Pi,lPi,3^i{l) 

r,(2)  +  Yi(  3) 

Pi^Pi^Ziit) 

%’)  +  ..•  + 

Yi(k)  +  ...  +  Yi(l)  pij . . .  ptik  . . .  pijzi(t) 

K(i)  +  y;(2) 

+  •  •  •  +  Yi(t)  piApi,2...phiZi(t) 

We  will  approximate  Ri{t) by  a  random  variable  of  Gaussian  (normal) 
type  with  the  same  mean  and  variance. 

Let  Yij  =  E[Yi(j)]  and  Ri(t)  =  E[Ri(t)\  where  E  denotes  expectation 
with  respect  to  5,(<)’s,  then 
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Ri(t)  =  Zi{t){Yi,ipi,i  +  >'i,2P«,2  +  •  •  •  +  Yi,tPi,t 

+  •  •  •  +  [Vij  +  . . .  4-  Yt,k  +  . . .  +  Yij]pij  . . .  piyk  ■  •  ■  Pi,i 
+  •  •  •  +  [K',1  +  •  •  -  +  •  •  •  Pi,«} 


and 

E[Rt(t)2} 


+  Vi#*  +  •  •  •  +  Y?tpi,t 

+  [y'j  + .  ■ .  +  Y*k  + ...  +  Y?t]pij . . .  j pi<k . . .  pi,i 
+  •  •  •  +  [VJj  +  •  •  •  +  Y?t]pi'i . . .  piyt} 


From  the  above  equations,  we  can  calculate  the  variance  of 
Var (R,(t))  =  E  [jfc(t)  -  Rjj)}2  =  E[Rx{t )2]  -  Rji)2 


The  inventory  available  to  fill  sales  for  period  t,  i.e.  the  starting 
inventory  for  period  t,  denoted  by  Ii(t)  can  be  written 


m  =/,(d  +  ft((  _  i)  -  g  s,(S) 

S=1 

Since  the  quantity  Jj(l)  is  a  constant  and  we  know  the  distributions 
for  R,(s)  and  5,(s)we  could  obtain  the  exact  distribution  of  I,(t)  from  the  con¬ 
volution  of  the  distribution  of  /?,(s)  and  —  5,(5).  Howerver,  since  we  make 

the  above  approximations,  we  need  only  work  with  the  means  and  variances. 


3.5.1  Constraints  for  Meeting  Requisitions  of  Important  Items 

The  important  items  that  we  mentioned  in  the  assumptions  section 
are  “MICl”  items.  MIC  is  an  acronym  for  Mission  Impact  Code.  If  a  shortage 
of  a  item  has  ever  caused  a  mission  to  be  cancelled,  the  item  is  coded  as 
‘“MICl.”  Other  items  are  coded  “MIC2”,  “MIC3”,  or  “MIC4.”  We  require 
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that  each  period’s  requisitions  (Sales)  for  MICl  items  be  met  with  at  least 
probability  a,i.e. 

P{Ii{t)  >  Si(t))  >  a, for  ie  MICl  for  t  =  1,...,AT 

Let  Di(t)  =  Y?k= 1  Si(k),  then  we  can  rewrite  the  constraint  for  the 
MICl  items, 

P{Ii(t)  >  Si(t)}  >  a 

as 

P{Ii(  1)  +  Ri(t  -  1)  -  £  Si(s)  >  Si(t)}  >  a 

j=i 

or 

P{Ii(  1)  +  Ri(t  -  1)  >  Di(t)}  >  a 

Taking  the  5,(<)’s  to  be  normally  distributed,  employing  linear 
stochastic  decision  rules  and  approximating  the  i?<(t),Si(t)combinations  by 
normally  distributed  variates  (e.g.  using  their  means  and  approximating  vari¬ 
ances),  we  will  invert  these  chance  constraints  in  the  process  of  obtaining  a 
convex  nonlinear  programming  problem  to  determine  the  coefficients  of  the 
optimal  linear  stochastic  decision  rules. 

3.5.2  Constraints  for  Matching  “due-ins”  to  “due-outs” 

We  call  an  unfilled  requisition  a  “due-out”  and  each  item  of  an  out¬ 
standing  order  a  “due-in.”  For  MICl  items  (high  priority  items),  we  wish  to 
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have  each  “due-out”  covered  by  “due-ins.”  So  for  these  items,  the  sum  of  all 
Sales  before  period  t  should  be  less  than  or  equal  to  the  initial  inventory  plus 
all  orders  through  the  current  period  t,  i.e. 

El=,  Si(k)  <  7,(1)  +  ELi YiW  for  t  =  1, . . . ,  N  . 

If  we  let  0{(t)  denote  the  cumulative  orders  through  period  t,  i.e. 
Oi(t)  =  i  (hence  Yi(k)  =  0{(k )  —  0{(k  —  1)),  then  we  wish 

to  have 


Di{t  —  1)  <  7.(1)  +0,(<) 

These  probability  one  chance  constraints  requirements  can  be  met  by 
constraints  on  the  stochastic  decision  rules  employed  for  the  ^-(fcj’s  as  induced 
by  these  constraints  on  the  0,(0’ s.  For  example,  if  these  are  to  be  linear, 

Oi(t)  =  ai(t)Di(t)-bi(t)li(  1) 

with  constraints  a,(f)  >  1  and  0  <  6,(2)  <  1  for  MICl  items  will  insure  proba¬ 
bility  one  satisfaction  of  these  MICl  constraints. 

3.5.3  Constraints  on  Dollars  Obligated  in  Each  Period 

Since  the  amount  of  money  permitted  to  be  invested  in  the  inventory 
is  limited  to  a  budget  amount,  ImaXi  we  must  try  to  insure  with  high  probability 
that  the  optimal  Yi(t)  values  lead  to  investments  in  inventory  that  are  less  than 
or  equal  to  the  budget  limit. 

* 

As  noted  earlier,  dollars  are  obligated  for  each  order  as  soon  as  the 
order  is  placed  and  funds  are  received  for  the  amount  of  each  sale  (requisition) 
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in  the  same  period  as  the  sale,  even  if  we  are  unable  to  fill  the  sale  from  stock 
on-hand  at  the  beginning  of  the  period. 

The  money  obligated  for  item  i  at  period  t  then  is 

c,{/,(l)  +0,(<)-A(*)} 

The  money  obligated  for  all  items  then  would  be 

+O.W-AW) 

i= 1 

So  our  dollars  obligated  chance  constraints  can  be  formulated  as 
P{£  [Ci(/<(1)  +  0,(t)  -  Di(t ))]  <  Imax}  >cj>  for  t  =  1, . . . ,  N 

t=i 


where  4>  is  close  to  one. 
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3.5.4  Summary  of  the  Chance  Constrained  Problem 

Using  informationally  feasible  linear  stochastic  decision  rules,  assum¬ 
ing  the  Si(t) are  independent  Gaussian  random  variables,  and  approximating 
the  i2,(t)by  Gaussian  random  variables  as  indicated  before,  we  will  have  the 
following  constraints  and  stochastic  expressions  to  invert  and/or  re-express  in 
terms  of  the  coefficients  of  the  stochastic  decision  rules.  We  then  schematically 
form  the  convex  programming  problem,  with  linear  and  quadratic  inequality 
constraints  in  terms  of  these  coefficients  and  an  objective  function.  Solution 
of  this  problem  will  determine  the  coefficients  in  the  linear  stochastic  decision 
rules  which  will  be  an  optimal  vector  of  linear  stochastic  decision  rules  for  this 
chance  constrained  programming  problem.  To  summarize  where  we  are  so  far 
in  formulae, 
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max  EEf(fi(  1)  +  Ri(t  - 1)  >  A(()} 

t=l  t=l 

s.t. 

P{Ii(l)  +  Ri(t  -  1)  >  Di(t)}  >  a  Vie  MICl; 

for  t  =  1, . . . , N 

P{YZ=i  [ci(/,(l)  +  Oi{t)  -  Di(t))]  <  Imax}  >  4>  for  <  =  1 , . . . ,  A’ 

Oi(t)Di(t)  -  bi{t)li(  1)  =  Oi(t )  for  i  =  1, . . .  ,n  ; 

for  *  =  1, . . . ,  N 
di(t)  >  1  Vie  MIC1  ; 

for  t  =  1 , . . . ,  Ar 
<ii(t)  >  0  for  i  ji  MICl  ; 

for  t  =  1 , . . . ,  N 

0  <  6,(0  <  1  Vic  MICl  ; 

for  t  =  1 , . . . ,  A’ 
0  <  6,(0  for  if  MICl: 

for  t  —  1 , . . . ,  Ar 

3.6  Inverting  the  Chance  Constraints 

Next,  we  invert  the  chance  constraints  to  form  an  equivalent  deter¬ 
ministic  convex  programming  problem. 

3.6.1  Inverting  the  Constraint  for  Meeting  Important  Items’  Req¬ 
uisitions 

We  can  rewrite: 


P{/,(1)  +Ri(t-  1)  >  D{(t)}  >  a 
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P{Ri(t  -  1)  -  Di(t)  <  7,(1)  }  >  a 


From  3.5,  we  know  the  means  and  variances  of  the  Gaussian  random  variables 
with  which  we  will  approximate  the  distributions  of  the  Ri(t  —  l)’s.  Further, 
since  we  have  the  mean  and  variance  of  the  S;(t)’s  from  their  joint  distri¬ 
bution  (assumed  known),  we  can  easily  calculate  the  mean,  £>,•(<),  and  vari¬ 
ance,  Var (Di(t)),  of  the  £),(<)’ s.  By  the  linearity  of  expectation,  the  mean  of 
Ri(t  —  1)  —  Di(t)  is  then 

Ri(t  -  1)  -  Di{t)  =  Ri(t  -  1)  -  D~(J) 

and  the  standard  deviation,  denoted  a ,(t),  is 

*i(t)  =  y/VarWt  -  1))  +  Var(D,(0) 

Since  the  random  variables  have  Gaussian  or  “normal”  distribution, 
the  resulting  distribution  is  Normal  with  mean  Ri(t  -  1)  -  D,(t)and  standard 
deviation  To  transform  into  a  Standard  Normal  expression,  we  subtract 

the  mean  and  divide  through  by  the  standard  deviation  in  the  chance  constraint 

n r Ri(t  - 1)  -  Di(t)  -  R,(t  -  1)  -  DM)  .  /,( 1)  -Ri(t-  Vj-njj), . 

— ^<0 - 5 - W) - 1 

Let  Ar(-)  denote  the  standard  normal  cumulative  distribution  func¬ 
tion,  then  our  constraint  becomes 
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But  A/”  (•)  is  a  strictly  one-to-one  function  and  therefore  has  an  inverse 
which  we  denote  A r~x  (•).  Also  note  that  a  >  .5  implies  that  Af-1  (a)  >  0. 
Applying  Af-1  (•)  to  both  sides,  we  get 

7,(1)  -Rt(t-l)-Djj)  ^  ir_w  x 

- W) - -A  (o) 

Since  <7,(t)must  be  positive,  this  implies 

7,(1)  -  -  Djfj  >  ai(t)Ar~1  (a)  >  0 

By  introducing  spacer  variables,  Vi(t),  we  have 

7,(1)  —  Ri(t  —  1)  —  D,(t)  >  v i(t)  >  <7,(<)A/’-1  (o)  >  0  and  hence 
7,(1)  —  Ri(t  —  1)  —  D{(t)  —  t\{t)  >  0  and  u,(<)  >  0  and  implying 

[v,-(t)]2  >  [<T,(<)]2[Ar_1  (a)]2-  Therefore  0  >  [<r,(t)]2[A^— 1  (a)]2  —  u,(<)2. 
(Note  (< <Ti(t )]2  =  Vu(Ri(t  -  1)  -  Di(t)).) 

Combining  these  results  with  the  results  for  the  mean  and  variance 
of  Ri(-)  from  3.5,  the  chance  constraints  for  MICl  items  is  inverted  to  the 
following  deterministic  constraints 

v,{t)  >  0 

7,(1)  -/?,(<-  1)  -  Dji)  -v<(t)  >  0 
(o)]2  -  Vi(t)2  <0 

Vie  MICl  ,for  t  =  1, . . . ,  A;  . 

where  <t,(<)2  is  Var {R,(t  —  1)  —  D,(t)). 
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3.6.2  Inverting  the  Constraints  for  Dollars  Obligated  in  Each  Period 
The  chance  constraints  for  dollars  obligated  each  period, 

+  Oi(t)  -  Di(t))]  <  Imax}  >  <f>  for  t  =  1 , . . . ,  jV 

tsl 

can  be  rewritten  as 

fMEMW)  )1  +  E[c.(Oi(<)  -  AM))  <  I  max}  >  <t> 

1  =  1  i=l 

then 


^{Eh(0.(<)  -  £>.(!))]  <  Imax  -  Elct/itl)  )]}  >  * 

1  =  1  1=1 

Using  our  stochastic  decision  rule  restriction  0,(f )  =  ai(t)Di(t)  —  &,(<)/,(  1)  the 
constraint  becomes, 


e{EM<*(0A(0-M0/i(i) 

i=l 


-  A(0)1  <  Imax  ~  Efet^l')  )D  5  * 
1  =  1 


or 

e{E[c,(Z),(i)(°.«)  -  1)  -  mh(  1)  ))  <  Imax  ~  E[C,(/,(1)  )]}  >  <> 

1  =  1  1  =  1 

or 


e<EMA(<)M0  - 1))]  <  Imax 

1=1 


+  EKM<)  -  Cij/.d)  ]} 

1=1 


><t> 


By 

M 

the  linearity  of  expectation,  the  mean  of  H"=i{c,(Z?,(f)(a,(<)  —  1)]},  denoted 
will  be 
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MO  =  £?=i {ci[(ai<  -  l]Di(t)} 

where  Dt(0  denotes  the  mean  of  the  cumulative  sales  which  we  can 
obtain  from  the  known  joint  distribution  of  sales. 

Similarly  by  stochastic  independence  of  Di(t)  for  different  i, 

Var(E,= ,  -  1)]})  =  EUMMO  -  l)JVar(A(i))] 

Let  <5(0  represent  the  corresponding  standard  deviation.  Using  the 
Gaussian  approximation  to  the  D{(t)  distribution  as  before  so  that  we  can 
reduce  to  standard  normal  expressions,  we  can  rewrite  this  chance  constraint 
as 

/{ - m - - - m - 1  £ 

Letting  A'  (•)  denote  the  standard  normal  distribution,  gives 


Using  the  same  argument  as  in  3.6.1  with  <j>  >  .5,  we  then  invert  this 
constraint  to 


/max  +  O6‘(0  -  C,)/,(l)  ]  -  fl(t)  >  (<f>)  >  0 

i=l 


Introducing  spacer  variables,  w(t),  we  get 


/mar  +  J^KMO  ~  C,)/,(l)  ]  -  n{t)  >  w(t)  >  S( t)Ar~l  {(f>)  >  0 

i=l 


•  •  yielding  the  following  equivalent  deterministic  constraints, 
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w(t )  >  0 

/max  +  E?=1[(W-C, )/<(!)]  -MO  "WO  >  0  (3.2) 

PftlAT-1  {<j>) ]2  -  iy(t)2  <  0 

for  <  =  1, . . . ,  N  where  /z(<)and  <52(t Represent  the  mean  and  variance 

ofE?=i[c,(A(OMO-i))]- 


3.6.3  Re-expressing  the  Functional 


The  objective  of  the  chance  constrainted  program,  as  stated  before, 
is 

N  n 

m*x EE />{/<(!)  +  Ri(t  -  1)  >  A(0) 

t=l  t=l 

We  would  like  to  maximize  this  sum  of  probabilities,  but  doing  this  would 
lead  to  a  very  complicated  problem  for  which  no  convenient  or  constructive 
solution  method  is  available,  so  instead  we  propose  the  following  two  surrogate 
expressions  whose  maximumizations  are  coherent  with  maximization  of  these 
probabilities. 

(1)  Maximize  the  minimum  of  these  probabilities. 

(2)  Maximize  the  sum  of  the  expected  values  of  the  random  variables, 
i.e.  max  X)f=iE”=i [/»(!)  +  Ri(t  —  1)  —  A(0]  subject  to  conditions  implying 
that  these  terms  are  non-negative. 


Maximizing  the  minimum  of  these  probabilities.  If  we  choose  to  maxi¬ 
mize  the  minimum  of  these  probabilities,  then  in  an  exactly  analogous  manner 
to  section  3.6.1,  we  rewrite 
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P{Ii(  1)  +  Ri(t  -  1)  >  Diit)} 


as 


7?,(*  -  1)  -  D,(t)  -  Rj(t  - 1)  -  A(0  7,(1)  Dj(t) 

<7i{t)  ~  ax{t) 


to  get  to  the  standard  normal  distribution. 

Let  7,(t)be  defined  such  that 

M  =  JV (7,(0) 

so  MibM^m=7l(t, 

Since  Ar  (•)  is  a  strictly  increasing  function,  maximizing  on  the 
7,(t)will  maximize  on  the  A^(7 i(t))  and  hence  on  the  probabilities  of  filling 
the  requisitions  when  placed.  Maximizing  on  the  minimum  of  these  probabili¬ 
ties  can  be  stated  as 


max  7 

s.t.  (3.3) 

7 i(t)  >  7  for  i  =  1, . . .  ,n  ,for  t  =  1, . . .  ,N 


or 


max  7 
s.t. 


>  7  ,ori=i . „ 

for  t  =  1 , . . . ,  Ar 


(3.4) 
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We  re-express  this  new  constraint  in  terms  of  ai(-),bi(-),and~i‘  as  fol¬ 
lows. 

Note  0  <  7  <  1  since  7  represents  a  probability.  So, 

7,(1)  -Rj{t  —  1)  —  Dj(t)  >  .  q 

<n(t)  ~  ~ 

Since  <7,(t)  >  0,  we  can  rewrite  this  expression  as 

7,(1)  -  Ri{t  -  1)  -  DJfj  >  ffi(t) 7  >  0. 

Introducing  spacer  variables,  U{(t),  we  have 
7,(1)  -  Ri(t  -  1)  -  Di(t)  >  Ui(t )  >  cri(t) 7  >  0. 

Yielding  the  following  deterministic  constraints, 

Ui(t)  >  0 

7,(1)  —  Ri(t  —  1)  —  A(7j  —  Ui(t)  >0 

^,(0272  ~  Ui(t)2  <  0 

for  i  jt  MICl  ,  for  t  =  1, . . . ,  Ar  . 

where  <7,(t)2  is  Var(7?,(f  —  1)  —  A(0)- 

Maximizing  the  sum  of  the  expected  values  of  the  original  objec¬ 
tive  function.  If  we  choose  the  surrogate  of  maximizing  the  sum  of  the  ex¬ 
pected  values  of  the  random  variables,  then  our  new  objective  function  would  be 
maxHtliirr=i[-^(l)  +•#,(*  —  1)  —  Z>,-(<)]  subject  to  the  following  additional  con¬ 
straints  to  insure  each  these  terms  are  non-negative  7,(1)  +  Ri[t  —  1)  —  D,(t)  > 
0  for  i  Ji  MICl 
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3.7  Summary  of  Equivalent  Deterministic  Convex  Pro¬ 
gramming  Problems 

We  now  combine  the  results  from  section  3.6  to  obtain  the  following 
convex  deterministic  non-linear  programs.  All  the  definitions  are  as  before. 

3.7.1  Using  Surrogate  Objective  (1) 

If  we  choose  surrogate  objective  (1),  to  maximize  the  minimum  of  the 
probabilities,  we  obtain  the  following  convex  deterministic  program. 
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max  7 
s.t. 

Ui(t)  >  0  for  i  /  MICl 

for  t  =  1 , . . . ,  Ar 

/,( 1)  -  Ri{t  -  l)  -Di{t)  -  Ui{t)  >  0  for  *7  MICl 

for  t  =  1 , . . . ,  N 

o,(0V  <  0  for  i  /t  MICl 

for  t  =  1 , . . . ,  N 

7,(1)  -  f?<(<  -  1)  -  ~Di(t)  -  Vi(t )  >  0  Vie  MICl  ; 

for  t  =  1, . . . ,  Ar 

Var(7Z<(t  -  1)  -  A(<))[.A7_1  (a)]2  -  t>,-(<)2  <  0  Vie  MICl  ; 

for  t  =  1 , . . . , 

Vi(t)  >  0  Vie  MICl  ; 

for  t  =  1, . . . ,  N 

-f  L?=1[(b,(t)  -  c,j/,(l)  }  -  fi(t)  >  0  for  t  =  1,...  ,N 

S2^)^-1  (4>)]2  -  w(t)2  <  0  for  t  = 

w(t)  >  0  for  /  =  1, . . . ,  N 

ai(t)  >  1  Vie  MICl  ; 

for  t  =  1 , . . . ,  N 
ai{t)  >  0  for  i 4  MICl  ; 

for  t  =  1 , . . . , N 
0  <  6,(0  <  1  Vie  MICl  ; 

for  t  =  1 , . . . ,  N 
0  <  6,(0  for  i  /  MICl  ; 

for  t  =  1 , . . . ,  Ar 
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where  62(t)is  the  variance  of  £”=1 [c,-(JDj(/)(a,-(<)  —  1))]  and  [<7,(t)]2  is 
the  variance  of  Ri(t  —  1)  —  Di(i). 

3.7.2  If  we  Choose  the  Surrogate  Objective  (2) 

If  we  choose  surrogate  objective  (2),  maximizing  the  sum  of  the  ex¬ 
pected  values,  then  we  have  the  following  convex  deterministic  program. 
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max5m  i*0-)  +  Ri(t  - 1)  -  Di{t) 

t=l  1=1 

s.t. 


7,(1)  +Ri{t-1)-Di(t) 

> 

0 

for  i  £  MIC1  ; 
for  t  =  1 , . . . ,  N 

7,(1)  -Jfc(t-l)-  Di(t)-Vi(t) 

> 

0 

Vie  MIC1  ; 

for  t  =  1 , . . . ,  N 

Var (Ri(t  -  1)  -  Di(i))[M-'  (a)]2  -  Vi(t)2 

< 

0 

Vie  MICl  ;  • 

for  t  =  1, . . . ,  AT 

Vi{t) 

> 

0 

Vie  MICl  ; 

for  t  =  1, . . . ,  N 

imax  +  U=il(Ut)-Ci)ii(i)]-n(t) 

> 

0 

for  t  =  1 , . . . ,  N 

6\t)[Ar-'  (<0]2  -  w{tf 

< 

0 

i—H 

II 

-to 

V-» 

£ 

w(t) 

> 

0 

r-H 

II 

Vh 

<8 

0,(0 

> 

1 

Vie  MICl  ; 

for  t  =  1 , . . . ,  N 

Oi(0 

> 

0 

for  i/L  MICl  ; 

for  i  =  1, . . . ,  N 

0  <  bi{t) 

< 

1 

Vie  MICl  ; 

for  t  =  1 , . .  • ,  N 

0 

< 

m 

for  i  t  MICl  ; 
for  t  —  1 , . . . ,  N 

where  <52(<)is  the  variance  of  X3r=i[c<(A(0(G«(0“l))l  and  [°‘i(0]2  is  the  variance 
of  Ri(t  -  1)  -  Di(t). 
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3.7.3  Conclusion 

For  either  surrogate  objective  function,  a  convex  quadratic  program¬ 
ming  problem  is  obtained  for  <Zj(-),  &;(•),  7,  and  the  spacer  variables  in  these 
expressions.  The  optimal  solutions  for  these  will  lead  to  expressions  for  the 
optimal  stochastic  decisions  in  terms  of  the  Yi(t)  . 


BIBLIOGRAPHY 


[1]  Brooke,  A.,  Kendrick  D.,  and  Meeraus,  A.  (1988).  GAMS:  A  User’s  Guide. 
The  Scientific  Press;  USA. 

[2]  Charnes,  A.  and  W.W.  Cooper,  “Chance-Constrained  Programming,” 
Management  Science ,  Vol.  6,  No.  1,  October,  1959,  pp.  73-79. 

[3]  Charnes,  A.  and  W.W.  Cooper,  Management  Models  and  Industrial  Ap¬ 
plications  of  Linear  Programming ,  Vols.  I  and  II,  John  Wiley  &  Sons,  New 
York,  1961. 

[4]  Charnes,  A.  and  W.W.  Cooper,  “Deterministic  Equivalents  for  Optimizing 
and  Satisficing  Under  Chance  Constraints,”  Operations  Research ,  Vol.  11, 
No.  1,  pp.  18-39. 

[5]  Charnes,  A.,  W.W.  Cooper,  S.  Duffuaa,  and  M.  Kress,  “Complexity  and 
Computability  of  Solutions  to  Linear  Programming  -Systems,”  Interna¬ 
tional  Journal  of  Computer  and  Information  Sciences,  Vol.  9,  No.  6,  1980, 
pp.  483-506. 

[6]  Charnes,  A.,  W.W.  Cooper,  and  M.J.L.  Kirby,  “Chance-Constrained  Pro¬ 
gramming:  An  Extension  of  Statistical  Method,”  Optimizing  Methods  in 
Statistics,  Academic  Press,  Inc.,  New  York  and  London,  1971,  pp.  391-402. 

[7]  Charnes,  A.,  and  M.J.L.  Kirby,  “Some  Special  P-Models  in  Chance- 
constrained  Programming,”  Management  Science ,  Vol.  15,  1968,  pp.  33-48. 


35 


36 


K 

* 


[8]  Hadley,  G.,  and  T.M.  Whitin,  Analysis  of  Inventory  Systems.  Prentice- 
Hall,  Englewood  Cliff,  1963. 

[9]  Yielding,  G.  Interim  Report:  Demand  Forecasting  and  the  Aggregate  Model 
Air  Force  Logistics  Management  Center,  Gunter  AFB,  AL  (in  print). 


